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SUMMARY 


original page is 

OF POOR QUALITY 


The radar response from the Amazon rain forest has been studied to 


determine the suitability of this region for use as a standard target 


to calibrate a scatterometer like that proposed for the National Oceanic 


Satellite System (NOSS). Backscattering observations made by the 


SEASAT-1 Scatterometer System (SASS) show the Amazon rain forest to be 


a homogeneous, azimuthally-Isotropic, radar target which Is Insensitive 


to polarization. The variation with angle of incidence may be adequately 


modeled as “ a0 + b with typical values for the incidence-angle 

coefficient, a, from 0.07 - 0.15 dB/deg, Variation of the intercept, b, 


was from -0.6 to =4. A dB. A small diurnal effect occurs, with measure- 


ments at sunrise being 0.5 dB - 1 dB higher than the rest of the day. 


Max i mum- 1 ike 1 i hood estimation algorithms presented here permit deter- 


mination of relative bias and true pointing angle for each beam. Specific 


implementation of these algorithms for the proposed NOSS scatterometer 


system is also discussed. 


Much of the information in this report can be found in the paper by 
I.J. Birrer, E.M. Bracalente, G.J. Dome, J. Sweet and G. Berthold, "o° 
Signature of the Amazon Rain Forest Obtained from the Seasat Scatterometer," 
IEEE Transactions on Geoscience and Remote Sensing , vol . GE-20, no. 1 , 
January 1982 , pp. 11-17. 
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1.0 INTRODUCTION * 

The National Oceanic Satellite System Scatterometer (NOSS SCATT) was 
designed to determine wind vectors on the ocean's surface by measuring 
radar backscatter from the surface. The SCATT was> to be a fan~beam scatter - 
ometer utilizing six dual "polar I zed (V and H) antennas. 

To achieve maximum wind-vector accuracy, careful cross-calibration 
between antennas must be performed. This report discusses the development 
of a calibration technique utilizing measurements of the scattering coeffi- 
cient a 0 from the Amazon rain forest. A brief discussion Is presented of 
the method used to discriminate rain-forest signals from others. Maps (used 
to screen Brazil data) generated for this purpose are on a 0.25° x 0.25° 
grid for locating rivers and a 0.5° x 0.5° grid for identifying vegetation 
types,. 

The suitability of the Amazon rain forest for use as a calibration 
target is shown in Section 4.0. From July to October 1978 numerous back- 
scatter measurements were made by the SEASAT-1 Scatterometer System (SASS) . 
These measurements show, the rain forest to be a? homogeneous, azimuthal lyr 
isotropic radar target which is- insensitive to polarization. Variation of 
0 ° with respect to incidence ang le may- be adequately, modeled' with a straight 
line fit of the form: , R \, a0 + b , where 0 is the? angle of Incidence 


(deg.). Measured values for the? incidence angle- coefficient, ay, range from 
0.07 _ 0.15 dB/deg depending on. beam and time of day. Values, for the 
intercept, b n range from- -O'. 6 to -4, 4- dB. A smal l- di urnaT effect was 
observed with measurements, made at sunrise- (0500-0630) be i ng' consistentl y- ,f 
0.5 - 1.0 dB higher than the rest of the day* 

To cross-cal ibrate the beams on- a mul tibeam scatterometer, twoi para** 
meters, a relative bias, and the true- aptqnna. pointing angle, must be deter- 
mined for each beam. Development of two maximum - 1 i kel ihood-est imatf on 
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algorithms to determine these parameters using Amazon rain forest data 
is discussed in Section 15.0. The first algorithm Is used to determine 
both relative bias and true pointing angle. A second and considerably 
simpler version Is then developed to monitor long**term changes In rela - 
tive bias due to drift In the Instrument gain and transmitter power. 
Specific processing required to Implement these algorithms for the SCATT 
Is also discussed in Section 5.0. 


-:?1 


2.0 BACKGROUND 

During the SEASAT program, the need for In-f light cross-cal ibration for 
antennas from a multibeam scatterometer became apparent. The originally 
proposed procedure using aircraft underf lights proved Inadequate. Based on 
the consistency of o° data over the rain forest obtained during SKYLAB, 

R.K. Moore of The University of Kansas suggested using the Amazon rain for- 
est data as a means of calibration. 

The accompanying map (Figure I) shows the vast extent of the Amazon 

2 

rain forest. It Includes an area of more than 3,000,000 km . Because of 
its equatorial location, seasonal effects are minimized and vegetation Is 
always present. The region is very flat and uniformly forested except for 
breaks In the canopy over the Amazon River and its principal tributaries. 

The predominant climate is hot and humid, with annual precipitation as high 
as 3000 mm. A slight decrease In rainfall is observed in the large central 
part of the rain forest during the months of October to December. A similar 
decrease takes place in the northern part during the months of June to 
September, in the northwest a more significant decrease in rainfall occurs 
during May and June with the lowest rainfall in July and August. 
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During the SKYLAB SL-2 mission vertical-polarization measurements were 
made over the Amazon rain forest at Incidence angles of 32° to 36°* Figure 


2 shows a sample time history taken during one of the SKYLAB passes over the 


eastern part of the forest. The scattering coefficient was stable with a 


maximum deviation of about 0.2 dB, For the combined set of SKYLAB measure- 


ments the mean value of a was reported as * -5.9 dB with a standard devla* 


tlon of 1 dB [1,2], This observed consistency prompted the further Invest I ■ 


gat Ion using the SEASAT scatterometer. 


3.0 DEVELOPMENT OF MAPS USED TO CLASSIFY BRAZIL TARGETS 


3. 1 Land-Water Map 


For the early analysis of SEASAT observations of Brazil, discrimination 


was achieved by using vegetation maps developed during the SKYLAB program [2], 


The 1° x 1° grid available was thought to be adequate. However, as noted 


in Section 4.2, the presence tn the SASS footprint of the Amazon River or 


one of Its principal tributaries markedly lowered the radar cross-section. 


Since the SCATT was planned to operate with an even smaller footprint, a 


finer map grid was needed to screen out data contaminated by the presence of 


rivers. 


For NOSS a digitized map was created with a 0.25° x 0.25° grid. The 


source map was produced by the Department of Cartography of the Brazilian 
Institute of Geography in 1971. This map of the entire Amazon basin Is 
scaled at approximately 3 Inches per degree. The digitized map was created 
by overlaying grids to subdivide the source map into 0.25° boxes. Codes 


were then assigned to each box using the following conventions: 
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FIGURE 2: Sample SKYEAB Time History over Amazon Rain Forest 
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0 * flatland only (no significant rivers) 

1 - some small rivers but mostly land 

2 - large rivers through box 

3 - rough terrain 

The results of this procedure were then encoded Into a computer algorithm# 

For a given latitude and longitude the algorithm will fetch the code of the 

box containing the location. This map was tested during the extensive analy- 
sis of SEASAT data discussed In Section h.O. 

For NQSS It was envisioned to use the computed location of the center 
of the Doppler cell for classification. A more careful screening could be 
performed by using the latitudes and longitudes of the corner points for 
a particular Doppler cell. Testing to be sure that all four corners showed 
a "land-only" condition would be a check against overlapping part of a 
river-filled box. 

3.2 Vegetation Hap 

To further classify Brazilian data a vegetation map for the Amazon basin 
was produced with a 0.5° x 0.5° grid. Codes were assigned to each box accord- 
ing to the conventions listed In Table 1. Six distinct types of vegetation 
are found In the Amazon basin. Types 1 and 2 are the forests which are suit- 
able for use as a standard target. Combinations of types were only designated 
If the dominant vegetation Was either; type 1 or type 2. Otherwise, they 
were coded with a 3# *** 5» or 6, depending on the dominant type. 
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This map was Incorporated into an algorithm similar to the land-water 
algorithm. For a particular latitude and longitude the code Is retrieved 
for the box containing that location. This algorithm was used during the 
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Type (Code) 

TABLE 1 

Codes Used In Vegetation Ma/> 
Description 

•0 

Outside of Brazil Ra la Forest 

1 

Humid Upper Amazon 

2 

Humid Firm Land 

3 

Humid Floodlands 

k 

Humid Plains 

5 

Mountainous Region 

6 

Savannah 

7 

Combination Types 1 ,k 

8 

2 ,k 

9 

1 > 2 , k 

10 

1,2 

11 

2,6 

12 

2 ,k ,6 

13 

1 , 2 , 3, 6 

Ik 

1,6 

15 

2 , 3,6 

16 

2,3 

17 

1,2,3 

18 

1,3 

19 

2,5 

20 

2 , 5,6 


SASS analysis to remove all data from vegetation types 3,k,5, and 6 . A 
similar approach Is recommended for a NOSS-type of system. 
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k . 1 The SEASAT-A Scatterometer System (SASS) 

The SASS operated using four dual-polarized (V and H) 0.5° X 25° fan- 
beam antennas. The antennas were pointed k5° from the satellite subtrack 
to give a "star-like 4 illumination pattern on the surface [ 3 ]. (see Figure 3) « 
This configuration provided for pairs of nearly orthogonal measurements of* * 
the earth's surface, The SASS instrument transmitted a 100-.wa.tt lk. 6 -GHz 
signal. The backscattered signal was spread out in frequency due to the 
Doppler effect. Using 15 parallel channels (each consisting of a Doppler 
filter, a square- law detector, and a gated integrator), the reflected power 
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The mean values of the 61 Integrated voltage levels of signal plus noise t 

i 

and noise alone from each of the 15 channels were transmitted to the ground. j 

■ ! 

The power reflected from the earth's surface was determined from these vol- 

! 

tage pairs, the cal ibrAted gain of the receiver, and the known Integration “ { 

time, as Indicated In Figure **. 

Figure 5 shows an enlarged view of one of the 15 Doppler cells synthe- 
sized along the beam. As Illustrated, the Doppler-cell area Is determined 
by the 0.5° 3-dB beamwidth in the narrow-beam plane and the Doppler filter 

\ 

along the beam. Due to the 1,89-sec measurement time the Doppler ceil Is | 

smeared in the alongtrack direction to give the resultant cell shown In the 
lower Inset. The integrated Doppler cells are between 16-20 km wide and I 

50-70 km long. j 

The frequencies selected for the Doppler filters result in three near- 
nadir resolution cells (0° - 13° incidence angle) and twelve off-nadir cells 
(20° - 65 ° incidence angle). Only the off-nadir measurements were used in 
this analysis since the near-nadir returns were outside the main antenna beam. 

The value of cr° can be determined by solving the radar equation in the 


following form: 


e -/W 


r 2 

<£•) . 

b 0 q°h(f) 


R m 3 


•R 




dA 


where: 


P-. = received power 
K 

P_ = transmitted power 
T 

X = transmitted wavelength 


L- = miscellaneous losses 



( 1 ) 




FIGURE k: Calculation of SASS Received Power 
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G 0 “ peak antenna gain 

■ relative gain at center of cell 

G 0 

h(f) * frequency response of Doppler filter 
ft ■ range to cell center 


Equation (l) may be solved for a 0 as: 

„ (W 3 p / 

0 - ■ ■¥J-T75tr- 

P r A L S fa 0 'G 0 ' A 


assuming: 


( 1 f o < f < f u 
h<f) - 1 

0 otherwise 


A » L R* a , where L Is the distance on earth's surface from f^ to f y 


f - center frequency 

B « noise bandwidth (Doppler filter) 
n 


f oU f « upper and lower Doppler**f I Iter cutoff frequencies 

Xj *’ u 

<f>^ = narrow dimension antenna beamwldth (radians) 


r _ II 

V f ° T 


f » f + — 
u o 2 


A particular cr° Is calculated by evaluating equation (2) using measured 

c 

values for P R and P T , values for G Q and ■— obtained from a lookup table, 
and computed values for R and A [4] . 


h . 2 Target Stability of the Amazon Rain Forest 

Initially the SASS data were processed by scientists at NASA Langley 
Research Center to produce time histories of o° over Brazil [4]. Figures 


(-)'♦« 

s\\oy{ samp! e time histories of Beam I vertical-polarization measure- 
ments made during an orbit Identified as Rev 952* The radar-cross-sect Ion 
stability over the Amazon rain forest was similar to that observed during 
the SKYLAB mission* When the rain forest results are compared with the 
observations over the ocean* three striking differences are noted: (1) a 0 

Is consistently higher over the rain forest (except at the smallest Incidence 

angle), (2) o° Is much more stable over the rain forest, and (3) the varia- 
tion in a 0 with Incidence angle Is much less. As can be seen In both Figures 
6 and 7, SASS cells which Included part of the Amazon or one of Its major 
tributaries were easily Identified by the sudden 1-3 dB drop In a 0 . 

Since the goal of the research was to determine the Amazon rain forest's 
suitability as a standard target, tests were made for both temporal and 
regional stability. Comparing data taken at the same time of day over the ■ 
same flight line throughout the mission provided a useful test for seasonal 
effects. Figure 8 shows the results of one of these comparisons. Approxi- 
mately 10 measurements were averaged at each Incidence angle to estimate the 
mean a 0 . The deviation of the mean o° Is negligible near the peak antenna 
gain ( 8 1 w M°) except for the September 18 measurements at 1)0°. This dis- 
crepancy may be due to atmospheric attenuation- At the extremes of the inci- 
dence-angle range the dev I at Ion is less than 0*5 dB* However, due to the 
short lifetime of SEASAT (99 days, July-October) , further research is needed 
to determine if any seasonal effects exist. The values of the mean o° for 
a> particular time of day were compared over numerous f light* l ines, to test 
for regional effects. No consistent discrepancies were observed oyer the i . 
rain forest. As a result of these tests,, the entire SEASWf ita in-forest- data 
set was analyzed as a unit. 
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FIGURE 8: Target Stability During SEASAT Mission 
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^.3 Incidence Angle Dependence and Diurnal Variation 

The SASS made measurements over the Amazon rain forest twice dally. 
During the early part of the mission, the Brazil crossings were between 
0500 and O63O (local time at center of area) and between 1630 and 1830, 
During the last month of the SEASAT mission, the crossings took place 
between 0900 and 1200 and between 2100 and 2^00. For the subsequent analy- 
ses the data were binned by polarization, beam, Incidence angle and time 
of day. Mean values and standard deviations of a 0 were computed for each 
of these bins. 

Plots and regressions of mean or 0 versus Incidence angle were made 
for each beam to describe the small Incidence-angle dependence observed 
In the time histories. The sample sizes for estimation of the means varied 
approximately from 10 to 100 measurements. Figure 9 shows the results for 
Beam vertical polarization. The late-morning and nighttime data agree 
well. The results from Beams 1 and 2 Indicate that the late- afternoon 
measurements were also consistent with the late-morning and nighttime data. 
On the other hand, the early-morning measurements were from 0.5 to 1.0 dB 
higher than data at the other times of day. 

As can be’ clearly seen from Figure 9, the basic angular' trend is 


adequately modeled by the equation 


a (dB) = 30 + b (3)' 

where 0 is the incidence angle. Typical values for the incidence" angle 
coefficient a, as determined from regression analysis, ranged* from 0.077 
to 0.159 dB/deg depending on beam and time of day. 
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FIGURE 3: Mean a ys incidence Angle; Beam 4; V-polarization 
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Equation (3) may be rewritten as 


- 0/0 

_o * u m * 

0 (real ) Ke 


where : 


An 10 


K «* e 


10 

a An 10 


Table 2 gives a complete summary of the values of a, b, K, 0 Q determined by 

it 

regression using (3) • From the summary In Table 2 some obvious differences 
are observed between beams. The data appear to be most consistent In the 
Incidence angle range of 30° - 50°. The four SASS antennas had been cross- 


can ora tea Using tne eleven arazi i passes avaiiapie to mmsh Lang i ey scien- 
tists soon after the 5EASAT mission [A] . Good agreement was achieved In the 
mid-angular range. The remaining differences highlight the need for the more 
sophisticated cross-calibration algorithm (using all of the SEASAT passes) 
developed here. 

The value of 0° estimated from the regression fit agrees well with the 
reported SKY LAB results. The SASS, however, shows much less deviation abcut 
the mean. Bars haVe been drawn on Figure S to show the scatter in the measure- 
ments. At incidence angles near the nominal pointing angle of , the 
standard deviations ranged from OVl to 0.2 dB. At the extremes of the incH 
dence angle range the standard deviations ranged from 0.3 to 0.6 dB. This 

ft . 

Both unweighted regressions and regressions weighted by number of 
cases were tried. Weighted regression results are reported here because " 
they take into account the higher density of measurements in the angular 
region of greatest measurement accuracy* In most cases the differences 
in the "a" coefficients between the weighted and unweighted regress i6n 
were wi thin the errors of the estimated coefficients. 
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TABLE 2 

Regression Parameters for the 
Incidence Angle Response of SASS Data over Brazil 


Model : * a6 + b 


a 0 (real) ■ Ke °, where K, 0 O are determined from a, b* 


Beam Time of 
# day* 


• 0.129 

• 0.158 

■ 0.112 

• 0.096 


R 1 2 3 4 xI00* 


0,669 

0,764 

0.519 

0.428 


33.67 
27.41 

38.68 

45.30 


Data points 


- 0.077 

-0.089 

- 0.135 

” 0.155 


0.435 

0.391 

0.663 

0.774 


56.31 

49.00 

32,17 

28.09 


• 0.148 

•0.156 


0.874 

0.767 


29.30 

27.91 


•0.078 


0.366 55.61 


- 0.103 

- 0.112 

“0.132 


0.582 

0.510 


42.29 

38.65 


0.595 32.95 


Time Code: 


1 « 0500-0630 

2 = 0900-1200 

3 “ 1630-1830 

4 = 2100-2400 


No data avail able 


’Squared multiple correlation coefficient 




S 
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Increase In the standard deviation Is probably due to a decrease In measure- 
ment precision caused by weaker signals at angles well away from the peak 

antenna gain, and also may be caused by small beam-point lng-angle errors I 

* * $ 

that have more Importance where the antenna gain varies more rapidly with j 

angle. * j > 

The regression equations determined from the data were used to estimate 
values of cr° at 40° Incidence angle which are plotted In Figure 10 as a 

| 

function of time of day for Beam 2. The early morning return Is 0.6 - 0.9 
dB higher than at the other times of day. Table 3 summarizes these estimates 
for all of the beams. When this effect was first noted In the original 
analysis of the first II passes, all of the data had been screened using 
GOES Imagery to eliminate cases containing significant cloud cover. In the 
present analysis the data were net screened. However, since this was the 
dry season few of the measurements should have been corrupted. 

Initially It was believed that, because of the random orientation of 
the scatterers In the uniformly vegetated rain forest, the values of o° 
would be Insensitive to polarization. While most of the measurements were 
vertically polarized, enough horizontal polarization data were taken to 
allow a comparison. On Figure 11 ahe plots f both V-pol ahd H-pol data for 
Beam I. As expected the results appear to be Insensitive to polarization. 

No noticeable azimuthal -angle variation was observed In the rain- forest 
data. This Isotropy can be demonstrated by comparing the 0900 - 1200 V-pol 
data from Beam 2 and Beam 4. While the beams differ 180° In azimuth, the * I 

measurements show agreement within 0.1 - 0*2 dB. This may be seen by over-' * 

laying Figures 9 and 1 1 . 

\ 




Local Brazil Time 


FIGURE 10: Estimated a 0 Frew Regression Fit 
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TABLE 3 

Summary of Diurnal Effect at 40° Incidence Angle 
Vertical Polarization 


Time of Day 

ESTIMATED 

SCATTERING 

COEFFICIENT (dB) 

Average 

Beam 1 

Beam 2 

Beam 3 

Beam 4 

1 ur 

All Beams 

0500 *• 0630 

-6.91 

-6.70 

-6.34 

-6.47 

-6.60 

0900 - 1200 

-7.49 

- 7*62 

-7.39 

-7.41 

1 

Jr- 

09 

1630- 1830 

-7.33 

-7.40 


M tm "m m — 

-7.36 

2100 - 2400 

-7.53 

-7.31 

-7.49 

-7-54 

-7.47 


4.4 SEASAT Summary 

The data obtained from SEASAT-1 have shown the potential of the Amazon 
rain forest for use as a standard target* Since the diurnal effects appear 
to be limited to early morning, data from other times of day could all be 
used together to build a data base for cal Ibratlon. The small scatter (0,1 - 
0.6 HB) In the mean values of a 0 Indicates the remarkable stability of this 
region and hence Its usefulness as a standard target. A baste character- 
ization of the Incidence-angle dependence can be made using a simple straight- 
line model of 0°^ vs 9* 

The differences between the regression coefficients shown In Table 2 
for the different beams cannot be due to actual 0° differences because Beams 
1 and 2 covered the same area, as did Beams 3 and 4. Hence these differences 
must be In antenna patterns or in pointing-angle differences. This points 
out the need for use of such a "standard target" in future multibeam scat- 


terometers. 


- 26 - 


5.0 DEVELOPMENT OF OFF- NADIR ANTENNA BIAS CORRECTIONS FOR NOSS 


5. 1 Description of NOSS Scatterometer System 


The NOSS Scatterometer System (SCATT) design Is for a second generation 
fan-beam scatterometer using experience gained from SEASAT (Section 1 ) . 


Figure 12 shows the proposed SCATT ground-track coverage. Besides the four 
beams used In SASS, two antennas are pointed at 65° from the satellite 
orbit plane to aid In Wind-direction alias removal. All of the antennas are 
dual -polar Ized (V and H). The SCATT Instrument transmits a 100-watt 14.0- 
GHz signal. The backscattcred signal Is Doppler filtered Into 60 resolution 


cells with each cell approximately 10 km x 10 km. The off-nadlr bias 


correction algorithm Is to be applied to the outermost 50 Doppler cells of 


each beam. 


5.2 Development of Maximum Likelihood (ML) Algorithms 
to Estimate Off-Nadir Bias Parameters 


During SEASAT analysis, Brazilian rain-forest data were first used to 


correct for antenna biases In the off-nadlr cells. Estimates for relative 


bias and pointing angle were determined by comparing data from a limited 


number of passes. For NOSS a' formalized algorithm was desired to allow 
automation of the estimation of these parameters by using all available 
ScAtT rain-forest data. The next subsection (5-2.1) describes in detail 
the derivation of a maximum likelihood (ML) approach to estimate relative 


bias and pointing angle. Once pointing angle is determined, a ML technique 


may be Used periodically to re-estimate relative bias. This algorithm Is 


described in Section 5-2.2. 


This ML approach is an adaptation of the approach used In the SEASAT 
Wind Vector Algorithm developed by F.J. Wentz (Remote Sensing Systems, Inc.) 




1 
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FIGURE 12: Proposed NQSS Footprint 
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5 . 2 . 1 Relative Bias and Pointing Angle Estimation 

For an operational system the radar parameters, P T and in equations 
(l) and (2) are specified in the design and carefully measured prior to launch. 
The relative antenna gain, |p-, Is determined from pre-launch antenna pattern 
measurements and the specified pointing angle, 0 , Although Py is monitored 
during operation, only relative changes In P^. can be measured. The true 
absolute value of P^ depends on pre-launch calibrations, However, as was 
observed with the SASS, the actual values for P^., G Q and 0 p after launch can 
differ from the design values due to errors In the pre-launch calibrations [**] , 
Effects of these differences at each Doppler cell can be seen by writing 
equation (l) for both the designed and actual values of P^., G Q , 0 p and the 
corresponding 0° inferred from equation (2). For a given Doppler cell 
(incidence angle): 


V ‘/Sl'Jl a° D A 


P - D P D H 

R (inr) 3 7 

where the subscript D refers to the designed values. 

p t^ G C 2[G/G O (0 d )1 a A A 
p _ A 1 A° 0 A P A 

R 5^7 

where the subscript A refers to actual values. 


Equating the above expressions and solving for a° in terms of a° 


y ields 


_ 0 _ a[G/G 0 < e A p» 0 A 

° D ~ K"/ 


P T G 0 

where a = relative bias = ~ 

P G 2 
T 0 
D D u 








; y 
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By substituting a standard reference target, at the dell Incidence 

angle, for the value of cr°, the relationship In equation (5) can be used to 

estimate the relative bias, a, and the actual pointing angle 0 . Since the 

A k 

measurements, o^, are noisy, max 1 mum- 1 ike 1 i hood estimation Is used to 

determine a and 0 . 

A p 

Formulation of the max l mum- 1 i kel I hood estimation may be given as: a 

vector of noisy observations, z, exists which depends probabilistically on 
a parameter of interest, B. The M.L. estimate of B Is then given by the 
condition 


(z/B) 


maximum 


where: 


B = b ML (z) 


f is the joint conditional probability density function. 




v: 


b^ L (z) is the max I mum- 1 ike 1 1 hood estimate of B. 

This approach can be applied to estimate a and 0 . with the following 

K 


assumpt ions : 


(1) Radar measurements are independent and may be described 
as a Gaussian process 

(2) The noisy observations, z, are the parameter 0° 

(3) The parameter of interest, B, is given by the relation: 


a [G/G 0 (6 )r 

— «• o® - B(a 0 .) 

[G/G- (9 J] 2 S A P 

d k 


{k) G/G q ( 6^) is known from pre-launch measurements. 

Applying the above conditions, the conditional probability In expression (6) 
may be re-written as: 




j 5535333^6^9 
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whe re : 


a°/B 


Ha 0 .} £ - !,N}/B(ot,0 n )] - 

£° A p 

u -[a° D - B(a,e J] 2 /2(Aa°) 2 

5 e * A 
£-1 __________ 


( 7 ) 


N / — 

ri ■ II /2ir (Aajl) 2 
£=■1 

N' ■« number of measurements 
Aa° = standard deviation of the measurements. 

A A 

Since the purpose Is to select values of a,0 that maximize expression 

A p 

(7) j several simplifications can be made< The quantity, n, may be Ignored 

Since it U not a function of the parameters a, 9 * The quantity Acr° is assumed 

A 

to Be constant since any changes In the standard deviation of the backscatter 
measurement aS a function of ground track OVef the rain forest appear to be 
Smalls Lastly, the logarithm of the likelihood function may be maximized, 

A A 

s i nee the same values 0,0^ Will maximize both the likelihood function and the 
jog-likelihood function* Therefore, the condition to satisfy is: 


1/2 l (a° fi - 8)^ * max 


£=l 


where B = B (ot , 0 ) s 
A P 

Since the function is stored as a look-up taWe anti V*0ry 

difficult to express analytical ly, the log -I I keiihOod expression IH (8) is 
maximized In two steps. First a migrating Search Is performed to approxi- 
mately center a 3 x 3 matrix about the solution, then the maximum of a 
bivariate Interpolating polynomial fit to points on this matrix is found 

A A 

to obtain a * 0 p . 


u 
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The search Is done by defining a log-1 Ikel I hood matrix, 9 1 ^ j » as follows: 


9 i,j ■ ' ,/2 J l |o 1 ti ' e< vy ]: 


where; 


» a Q + lAa; 1 » -1,0,1 


0 jp ■ e op + -l i6 p'> J ■ "M,! 


Since the actual values for P^. and Gq should be reasonably close to 
the design specifications, Initially a is chosen to be 1.0. Similarly, 


the design-specified pointing angle, 0 , is chosen for the initial value 

D p 


of 0 O p. Reasonable values for Aa,A0 p appear to be 0.2 and 1.0 respectively. 


The values for g. . are then compared to find the maximum. If the center 
• > J 


point, gQ q, is not the maximum, a new matrix is formed as follows: 


ot . = ct + i Aa; i = -1,0,1 
i max 


e jp = e p + J 4e p > J • - ’ i°. * 

J K max K 


where: 


a ,0 are the a.,0.„ corresponding to the maximum of the previous 

" ax P m ax ' J p 


matrix. 


Searching continues until a matrix is formed which has the maximum at the 


center. 


Estimates for relative bias, a, and actual pointing angle, 0 p , are found 


by fitting a bi-variate, quadratic interpolating polynomial to the values 


of g. j. Values of a and 0 p which maximize the polynomial are then deter- 


mined analytically. The interpolating polynomial is of the form: 


g(a,0 p ) = a(3? + b{5 + cy 2 + dy + e0y + f 


where: 


= (u-cxq )/Aot 


<V 9 P >^ P 






M>uaaaaM^ac5agacaBa 


a * Li» - g 0 0 + 




9q -V 9q 

C > - n „ + 

2 J 0,0 2 


d .2o^ + %L 


• - »O,0 ■ »l,0 ‘ 9<M ■*®M 


f - 9, 


where g. . are the values of the likelihood matrix. 

The maximum is found by taking the partial derivatives of the inter- 
polating polynomial and setting them equal to zero. When the resultant 


equations are solved for a and 9 ? we get 

a a a + M(ed- - 2bc) 

0 *ac - e 2 


(Ha) 


0 = 0_ 


A0 (be - 2 ad) 


4ac - 


where a,b,c,d,e are defined' as above, 

A A A 

Final estimates for a and 0- p are made by averaging the values of and 

A 

0 obtained for each Doppler cell. 

P ML 

5.2.2 Long-Term lion Storing of Relative Bias 

Periodic estimates of the relative bias during the mission wi! r l provide 
a means to correct for long-term drists in transmitter power, if the pointing 

A 

angle is not expected to change once 0 p is determined, a simplified version of 
the algorithm in Section 5.2.1 may be used to estimate relative bias. 






HR 
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i 


In this case the parameter of Interest contains only one unknown, 
the relative bias a. The same procedure followed In Section 5*2.1 Is used 
to solve this one-parameter case with a likelihood vector, gj» replacing 


the matrix 9j j* The Vector Is defined as: 


N 


7 l a?] 2 ; I --1,0,1 


o l n 
4 £-1 


(12) 


where: 


a n ■ the NOSS SCATT measurement 
JT 

a? * the standard target 


oi| “ a Q + lAa 

a n * previously determined relative bias 


Vector 9 j Is then fitted with a quadratic Interpolating polynomial 
of the form: 


where: 

e 


g (a) «* r£ 4 + s3 + t 


(a - a Q ) 


(13) 


Aot 

9-1 * 9 | 

r - — " 9 0 + T 

9_, 9, 

s — r + x 

t - 9 0 

gj ■ the values of the vector calculated In equation (12). 

By differentiating equation (13) and setting the result equal to zero, 

A 

the maximum likelihood value of a Is found as: 



4 


■ A i.' ' ", 


, ,»1 ■ 9-1 

“o - to( g., i!i - «(,"+ i? 


A final estimate of o may be determined by averaging the estimates of a. 


obtained for each Doppler cell. 


5.3 Implementation of Off-Nadir Bias Algorithm for NOSS SCATT 


The flow chart In Figure 13 shows the relationship between the off- 


nadlr bias algorithm and SCATT data processing. Earth-located o° data 


are screened to locate data over the Amazon rain forest. Then the Amazon 


map described In Section 3.0 Is used to remove measurements that Include 


either the Amazon river or one of Its principal tributaries. Remaining 


rain-forest data are then supplied to the off-nadlr bias algorithm. 


Actual estimates of the relative bias and pointing angle are performed 


off-line from the main SCATT processing. This Is necessary for two reasons: 


first, Brazil passes for several days are needed to create a sufficient 




data base for estimation; second, considerable human Interaction and evalu- 


ation Is needed to Interpret the results. 


Figure 14 shows a flow chart of the processing required in the off- 


nadlr bias algorithm, inputs to the algorithms are the SCATT 0 data from 


the rain forest and the corresponding brightness temperatures T 0 measured 
by the Large Antenna Multichannel Microwave Radiometer (LAMMR) . Co- located 


LAMMR data are used to flag and remove potential rain-corrupted SCATT 


measurements. A data base is then created for use with the maximum- 1 1 kel I* 


hood procedure described in Section 5*2. Estimates for relative bias and 


antenna pointing angle are computed for every combination of beam, ceil and 


polarization. Approximate error bounds are determined using the results of 


several estimations. The following subsections describe in more detail the 


processing required. 









FIGURE 1 



TO LEVEL II PROCESSING 


3: Proposed Processing of NOSS Scatterometer (SCATT) Data to 

Obtain 0° 






SCATT DATA 


LAMMR DATA' 



FIGURE 1 A: Processing Required for Off-Nadir Bias Correction 
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5*3.1 Co- Location of SCATT and L AMMR 

To use LAMMR 37 GHz data to flag rain-corrupted SCATT measurements, 

T 0 measurements must be collocated with Individual SCATT a 0 measurements* 
The Input data streams for both SCATT and LAMMR are assumed to be grldded 
Into 50 km blocks and assigned a set of coordinates (l,j). Co-location 
Is accomplished by: first, sorting both LAMMR Tg's and SCATT a 0, s Indepen- 
dently by grid coordinates; second, for each coordinate pair (1,J) contain- 
ing SCATT data the corresponding LAMMR grid point Is located. Within the 
50 km grid each SCATT cell Is matched with the closest LAMMR cell. 


5.3.2 


Ing Rain Cells 


Co-located LAMMR 37 GHz Tg's are used to flag potential ra!n-corrupt<id 
SCATT cells. The brightness temperature of the Brazilian rain forest Is 
very high due to high emlssivlty and surface temperature. However, In the 
presence of rain only the top portion of the rain cloud and not the 
forest Is seen. Hence It Is believed that the temperature measured should 
be lower when rain occurs than in clear sky conditions. Rain flags will be 
set whenever the measured temperature falls below a fixed cut-off temper- 

ature> W 

Figure 15 shows the form of the 37 GHz brightness temperature histogram 
expected from the Amazon rain forest. The lobe at the lower brightness 
temperature corresponds to measurements made In rain conditions, while 
the lobe at higher brightness temperature corresponds to clear-sky conditions. 
A preliminary value for T cUt may possible be made from SEASAT radiometer 
measurements. However, it will be Important to make a number of histo- 
grams with LAMMR data to confirm the assumed distribution shown in Figure 
15 and to properly set T 









FIGURE 15? Nature of Expected Amazon T R 01 strlbutlon 
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5.3*3 Data Base Creat ion 

Due to the large volume of data produced by the SCATT, estimation of 
relative bias and pointing angle Is performed by using mean o° values com- 
puted for each pass over Brazil. Individual clear-sky SCATT measurements 
are grouped by beam, cell number and polarization. At the end of each Brazil 
pass the mean value, standard deviation, and average Incidence angle Is 
computed for each group containing more than twenty measurements. These 
estimates are also binned by beam, cell number and polarization. Estimation 
of relative bias and pointing angle takes place when each mean value bln 
has data from at least 10 passes. Since the bins will fill at different 
rates depending on the location of the satellite subtrack, provision is made 
to store up to 20 values In each. 

5.3*4 Standard Target Creation 

The maximum likelihood procedure described In Section 5*2 requires a 
standard target, Og(0), defined over the entire off-nadlr Incidence angle 
range. Determination of Og will require consistent human evaluation and 
decision. 

Initially plots and regressions of the mean value of o° versus Incidence 
angle will be made from the data base described above. These plots and 
regressions will be compared with SEASAT results described In Section 4.0 
and Appendix B. If there Is good comparison In trend, the SEASAT model, 
a dB * a ® + will be adopted with an appropriate level shift. If the trend 
of the SCATT data Is significantly different, a new model determined from 
the analysis of several passes may be used. The final procedure used in 
creating the standard target will be determined once the Initial comparisons 





are made. 
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Early fn a future mission with the NOSS design will be the Implementation 
of the ML technique to determine both relative bias and pointing angle 
j described In Section 5.2.2, Once a suitable value for pointing angle is 

determined for each beam, the simplified algorithm described In Section 5-2.3 

H 

j; may be used to monitor long-range transmitter drift. Both of these algorithms 

ii 

j, can be implemented using the same Inputs so switching between algorithms 

j; should be trivial. 

| ( ft 

L In this implementation the noisy data, o n In equations (8) and (12) 

ij * 

are the mean values stored In the data base (Section 5.3.3). Separate 
estimates will be made for each combination of beam, cell and polarization. 
Repeated' estimates should be made using subsequent data. Final estimates 
and' bounds may be calculated by computing the means and standard deviations. 

6.0 CONCLUSIONS AND FUTURE RESEARCH 

Detailed analysis of SEASAT scatterometer data has confirmed the Suita'- 

? ' 

i b fifty of the Amazon rain forest for use as a standard calibration target. 

h 

fJ 

A Simple straight line model for incidence angle variation appears to be 1 
quite adequate. Useful approximations for the model parameters may be deter- 
mined from' the SASS data. The accuracy of these approximations Is limited by 
remaining biases between SASS antennas and the lack of independent raiii-fOres't 
measurements for absolute calibration. 

Maxi mum- likelihood estimation algorithms’ have been developed' 2 to correct- 

i 

If , . • * • 

| off-nadir ah’tenria'' bia'Ses : for spa'ceborne SCatte rbmeters by utilizing 2 rain- 

forest data. Applicat ion of' those al'gbr i thirls for the proposed' N0S§’ scatter- 
j; ometer was' straight-forward. Presumably these techniques could be adapted 

ii easi ly fop cal ibrati'ng other future systems. 


Three Important areas of research need to be undertaken to extend the 
usefulness of Amazon rain forest data to correct off-nadir antenna biases. 
First, measurements must be made throughout the year to determine the exis- 
tence of any seasonal effects. The SEASAT measurements were made during 
only the dry season. Second, research needs to be done to confirm the pro- 
posed algorithm to screen out data corrupted by thick clouds and rain. The 
multifrequency radiometer measurements made with the NIMBUS 7 satellite 
could potentially contribute valuable Information In both of these areas. 
Third, a calibration program using aircraft overflights Is needed to better 
determine the absolute value of 0 ° . 
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CBRZLND 
C 


C 

'C 

C 

•°0 


LAND/ WATER ROUTINE FOR BRZILIAN RAIN FOREST 

SUBROUTINE LDBRZL 
PURPOSE 

THIS ROUTINE LOADS THE LAND WATERMAP 


& 


INTEGER ARRAY (6,10,64), CODE 
REAL GLAT(6)/2. ,0. ,-2. ,-4. ,~6. ,-8./ 
REAL GLON( 1 0)/74* ,72, ,70. ,68. ,66. ,64. , 
62. ,60. ,58. ,56./ 

REAL LAT,LONG 


READ LAND/WATER MAP 


(2.-LAT)/2. ) + 1 
(74.-L0NG)/2. ) + 1 
(GLAT(ILAT) - LAT)/0.25 ) 
JTEP = IFIX( (GLON(JLON) - L0NG)/0.25) + 
KELE = ITEP*8 + JTEP 
CODE = 1 + ARRAY ( ILAT , JLON , KELE ) 


ILAT = IFIX 
JLON = IFIX 
ITEP = 


IFIX 


RETURN 
• END 


BRZVEG 


VEGETATION CODE ROUTINE FOR BRAZILAIN AREA 
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00000060 
00000070 
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00000090 
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00000110 
000001 20 
00000130 
000001 40 
000001 50 
000001 60 


00000460 

00000470 

00000480 

00000490 

00000500 

00000510 

00000520 

00000530 

00000540 

00000550 


1 



. DO 10 1=1,6 



00000170 


DO 10 J=1 .10 



00000180 


READ (01 ,900) ( ARRAY ( I , J , K) , K=1 ,64) 

00000190 

900 

FORMAT (4X, 818) 



00000200 | 

10 

CONTINUE 



00000210 





00000220 1; 


RETURN 



00000230 





00000240 


ENTRY LNDH20 ( TLAT , TLONC- , 

CODE) 

00000250 

1 

4 




00000260 

4 




00000270 

1 

4 

PURPOSE 



00000280 

1 

4 




00000290 | 

'S 

4 

THIS ROUTINE 

FETCHES 

A LAND/ WATER CODE 

00000300 

'I 

4 

CORRESPONDING 

TO AN 

INPUT LATITUDE AND 

00000310 

4 

LONGITUDE 



00000320 

"t 

4 




00000330 

/ 

INPUT ARGUMENTS 



00000340 

1 

/ 




00000350 

1 
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* 

LATITUDE (DEG) 
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EAST LONGITUDE (DEG) 
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OUTPUT ARGUMENT 
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CODE 
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SUBROUTINE LDVEG 


PURPOSE 


THIS ROUTINE LOADS THE VEGETATION MAP 
INTEGER VBGTBL(5,3,64),VGC0DE 
REAL LAT, LNG 

REAL VLAT(3)/2.0,-2.0,-6.0/ 

REAL VLNG(5)/74. 0,70. 0,66. 0,62. 0,58.0/ 

READ VEGETATION MAP 


r 


I 




ft 

I 


10 

900 

J 


jVEGCD 

1 

/ 


T 

/ 

1 

/ 

t 

/ 

i 

i 

i 

1 

/ 


DO 10 I = 1,5 
DO 10 J - 1.3 

READ ( 02 , 900 ) ( VEGT BL ( I , J , K ) , K= 1 , 64 ) 

CONTINUE 

FORM AT (4X, SI 5) 

RETURN 

COMPUTE CORRECT VEGETATION CODE 
ENTRY VEGCB( TLAT,TLNG, VGCODE ) 
PURPOSE 


THIS ROUTINE INPUTS A LATITUDE AND A 
LONGITUDE. WITH THESE A VEGETAI ON 
CODE IS DETERMINED 


r 

j 

j 

f 

r 


INPUT 


ARGUMENTS 


TLAT * 

TLONG * 


1 

t 

\ 

t 

f 

/ 

1 

t 

V 


OUTPUT ARGUMENT 
VGCODE 


LAT = TLAT 
LNG = TLNG 


LATITUDE (DEG) 

EAST LONGITUDE (DEG) 

VEGETATION CODE 


CONVERT TO WEST LONGITUDE 


LNG = 360.0 - LNG 


DETERMINE VEGETATION CODE' 


ILAT = IFIX( (2.0 - LAT) / 4.0) + 1 
I LNG = IFIX( (74.0 - LNG)' / 4.0) + t 

INDEXX = IFIX( (VLAT(ILAT) - LAT) / 0.5) 
INDEXY = IFIXC (VLNG(ILNG) - LNG) / 0.5) + 1 
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} 


INDEXZ * INDEXX * 8 + INDEXY 

VGCODE « VEGTBL(ILNG,ILA'D,INPEXZ) 

RETURN 

END 
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SUBMODULE 


IDENTIFICATION 

■ .»■ i . -> 

1 . 12.1 

1 . 12.2 


1.12.3 
1.12. A 
1.12.5 


1 . 12.6 


1.12.7 


1 . 12.8 


TITLE 

Co- location of SCATT with LAMMR Tg's Over the Amazon 

Rain Flags 

Data Base Creation 

Relative Target Creation 

Relative Antenna Gain Interpolation 

Standard a 0 Interpolation 


Relative Gain Table and Antenna Pointing Angle Corrections 
Long Range Monitoring of Relative Biases 
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1 . 12 . 1 ( 1 ) 
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SUBMODULE REFERENCE NO. 1.12.1 


SENSOR; SC ATT 


ORIGINATOR/PHONE: E.M. BRACALCNTE, I.J. BIRRER/80^-827“3631 , 913-86^-^836 

SUBMODULE TITLE: S C ATT- LAMM R 37 GHz Colocation over the Amazon 

SUBMODULE FUNCTION: This algorithm colocates the LAMMR brightness temperatures 

at 37 GHz with the SCATT data over the Amazon. 

INPUTS: 

(1) SCATT data over the Amazon (50 km data blocks) 

(2) LAMMR brightness temperatures at 37 GHz over the Amazon (50 km data blocks) 

OUTPUTS: SCATT data vector with colocated LAMMR Tg's appended to data 


vector. 


VIII. AUXILIARY DATA: TBD 


EXCEPTIONS: TBD 


PROCESSING REQUIRED: TBD 

COMMENTS: Since the calculation of brightness temperatures is a level II 

computation and this is a LEVEL I SCATT algorithm, provision must be made 
to feed back the LAMMR Tg's to this algorithm on a non-real time basis. 

The input data stream for both the SCATT and the LAMMR is assumed to be in 
50 km data blocks and will be assigned subgrid coordinates (i,j). 

For these coordinates the following procedure could be used for colocation: 
(a) Sort the LAMMR and SCATT independently by grid coordinates into two row 
ordered lists as follows: 


LAMMRLIST = T. T. T. ... T. T. T, . . . T, 
b , b , d b , b , b b 


1,1 1,2 1,3 


1,N 2,1 3,2 


SCATT LIST = SCATT, SCATT, .... SCATT, SCATT, SCATT ... SCATT 


1,1 1,2 


1 ,Q 2,1 2,2 


where T^ Is the 37 GHz brightness temperature of grid coordinate 


SCATT is the scatterometer data vector at grid coordinate 
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1 . 12 . 1 ( 2 ) 


(b) Colocate SCATT and LAMMR data by simply moving down the lists In 
parallel. For each SCATT append the corresponding T. . If more 

U ITJ 

than one brightness temperature exists, choose the closest. 



t 
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I. SUBMODULE REFERENCE NO: 1.12.2 

It. SENSOR: SCATT 

III. ORIGINATOR/PHONE: E.M. Bracalente, I . J. Bf rrer/804-827-3631 , 913-864-4836 

IV. SUBMODULE TITLE: RAIN FLAG 

V. SUBMODULE FUNCTION: Set a rain flag If the LAMMR brightness temperature at 

37 GHz Is below a critical value, T cu -, to Indicate the possibility of rain. 

VI. INPUTS: 

LAMMR T B 's at 37 GHz colocated with SCATT cells (Submodule 1.12.1) 

VII. OUTPUTS: 

Ra Inf lag (RFLAG) 

VIM. AUXILIARY DATA: TBD 

IX. EXCEPTIONS: TBD 

X. PROCESSING REQUIRED: 

(a) If T r , 7 Is less than T , then RFLAG = I (true), otherwise RFLAG = 0 (false). 

XI. COMMENTS: 37 

(a) This procedure may be designed to handle a cell at a time or all of 
the cells for a particular beam. 

(b) The value of T cut must be determined experimentally from LAMMR data. 

This could be done by plotting a histogram of LAMMR data over the Amazon 

taken on the first several passes. An Initial guess for T „ will be 

cut 

provided based on SeaSat-SMMR data. 
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M2. 3(1) 


SUBMODULE REFERENCE NO: 1.12.3 

SENSOR: SCATT 

ORIGINATOR/PHONE: E.M. Bracalente, ).J. Bl rrer/80it-827~363 1 , 913-864-4836 

SUBMODULE TITLE; DATA BASE CREATION 

SUBMODULE FUNCTIONS: This module creates a data base cf SCATT data over the 

Amazon rain forest. Data that Is not flagged for rain (Submodule M2. 2) Is 
sorted by beam, polarization, cell no. At the end of each pass over Brazil, 
the mean NRCS, standard deviation of NRCS, average incidence angle, and 
average antenna broadbeam angle for each combination of beam, cell, and 
polarization are computed. These results along with the number of points In 
each bln are stored In a semi-permanent fashion (tape, disk) for later pro- 
cessing. 

Plots of mean NRCS ( a °") versus cell no. are made for any combination of beam, 
polarization, etc. for which data exists. 


INPUTS: 


SCATT data over Amazon. 


Rain Flag. (Submodule 1.12.2) 


OUTPUTS: 


Key Parameters 


No. of Parameters 


Mean NRCS (7) 


Standard Deviation Aa c 


Count (N) 

Average Incidence Angle 9 
Average Broadbeam Antenna Angle e 
Pointer (P) 


Temporary data storage must be sufficient to store results from up to 20 passes 
or approximately 60,600 parameters. 


VIII. AUXILIARY DATA: None 


IX. EXCEPTIONS: If there are less than - 20 points within a bin, no processing 

is performed. 

X. PROCESSING REQUIRED: 

(a) Sort data by beam (i), cel 1 ( j ) , polarization(k) (see note 1). 
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(b) At the end of each pass, for each combination of l,J,k where there are 
more than 20 values, compute 


A - P 


i J,k 


. "| T k 

0 = __1 y’ J,K 0 

1 N | ^ j pgj'.j » J >k>n 


Aa f j,k,A 


U,k o 2 “o 

y (a 0 . . . ) z - N. . . o°. . 


N u,k - 1 


I >J >k,A " N. . 


IJ.k n-1 


LJ,k 

I e 


I > j *k,n 


j k 

1 v ,J,K 


IJ,k,A W,j >k n =i 


E * » j > k 


P. . , “ P. . . + 1 
r,j,k i,j,k 


where N. , . is the number of points In each i,j>k bln. 

(c) Plot a 0 versus cell number (j) for different combinations of beam (I) 
and polarization (k) 

COMMENTS: 

(1) Semi -permanent storage should be provided for all of the input SCATT ahd 
LAMMR data over the Amazon. 

(2) a 0 are assumed to be in ratio form. 

(3) Well-written computer code should only heed to store running sums during 
each pass and not all of the data values. 


Note 1: Other sorting parameters may be required. 
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SUBMODULE REFERENCE NO. 1.12./* 

SENSOR: SCATT 

ORIGINATOR/PHONE: E.M. Bracalente, I.J, Bi rrer/80i»-827“3631 , 913-864-4836 

SUBMODULE TITLE: Relative target creation 

SUBMODULE FUNCTION: This submodule takes the summary Amazon data (Means and 

Standard deviations) and creates a table of cr O (0) at 1° incidence angle 
steps for the off-nadir incidence angle data of the SCATT. This submodule 
will require considerable human interaction off-line. 

INPUTS: 

Mean, Standard deviations, average, incidence angles from submodule I . 1 2 . 3 • 

Plots of cr G versus cell no. from submodule 1.12.3. 

OUTPUT: 

Table of NRCS in one degree steps, over the off-nadir incidence.anc|1e range 
for both H and V polarizations. 

AUXILIARY DATA; TBD 
EXCEPTIONS: TBD 

PROCESSING REQUIRED: Since this task will require consistent human evaluation 

and decision, It will not be possible initially to fully automate this procedure. 
The following steps are envisioned at this point: 

(a) Compare data from each polarization and beam with SeaSat Brazil model 
(or other models of Brazil if available). This will be done by: 

(1) Examining the plots from. Submodule 1.12.3* 

(2) Regressing Mean NRCS (o°) against the general model cr° = F(6) 
determined from SeaSat. 

(b) If the results from (a) show a good fit with SeaSat model (except for a 
level shift), the SeaSat model will be adopted as the standard target. 

(Level shifted, if necessary). The beam which best fits the SeaSat model 
for each polarization will be used to determine the level shift. 

(c) If the data from the SCATT is best modeled by a different function than 
was used with the SeaSat, this model may be chosen to generate the standard 
table. 

COMMENTS: 

(1) The exact approach taken will probably be determined from a preliminary 
analysis of a few passes. 

(2) If there is good agreement with the SASS model or a satisfactory new model 
developed, it may be possible to automate the procedure. 


N.. 
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I. SUBMODULE REFERENCE NO: 1.12.5 

II. SENSOR: SCAT! 

III. ORIGINATOR/PHONE: E.M, Bracalente, I.J. Bf rrer/804-827“3631 , 913-86^-4836 

IV. SUBMODULE TITLE: Relative Antenna Gain Interpolation 

V. SUBMODULE FUNCTIONS: Use a three-point Interpolation to determine the 

relative antenna gain at the Input Incidence angle. 

VI. INPUTS: 

Mean Incidence angle (0)and mean antenna broadbeam antenna angle (F) corresponding 
to each 0 and value of pointing angle 0 

n p 

VII, OUTPUT: 

Relative antenna gain G/G q (0) 

VIII. AUXILIARY DATA: 

to table of relative gain values in Module 1.6 

IX. EXCEPTIONS: TBD 

X. PROCESSING REQUIRED: 

(a) e = F + (0 “ ® _) 

p n p 

(b) Find £ | , the Integer smaller than £ 

Z ■ 0 

(c) P « £ - e ( 

(d) G/G Q (Z) = fiflli G(e ,- 1 ) + ( 1 — P 2 ) G (e , ) +^i£±il G (e,+1 ) 

where G/Go(Z) Is the relative gain at an incidence angle of Z degrees. 

XI. COMMENTS: TBD 


/ • 
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I. SUBMODULE REFERENCE NO: 1.12.6 

II. SENSOR: SCATT 

III. ORIGINATOR/PHONE: E.M. Bracalente, I.J. Bi rrer/80V-827~3631 , 913-864-^836 

IV. TITLE: Standard CT° Interpolation 

V. SUBMODULE FUNCTIONS: Use a three-point Interpolation to determine the 

, standard NRCS of the Amazon rain forest at the input incidence angle. 

| i VI. INPUT: incidence angle (0) 

■- ■ VII. OUTPUT: Standard NRCS a° 

' i! VIII. AUXILIARY DATA: 

Tabular form of standard radar target 
or 

Tabular form of relative target from Submodule i . 1 2 . 4 . 

IX. EXCEPTIONS: TBD 

X. PROCESSING REQUIRED: 

(a) Find 0| the largest integer smaller than 0 

(b) P = 6 - ©j 

(c) a°(9) = or°(0, - 1) + (1 - P 2 )a°(0 | ) + £i£tiia°(e j + 1) 

where a°(Z) Is the NRCS of the Amazon rain forest at an incidence angle 
of Z degrees. 

XI. COMMENTS: TBD 
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SUBMODULE REFERENCE; 1.12,7 
SENSOR; SCATT 

ORIGINATOR/PHONE; E.M. Bracalente, I.J. Bl rrer/804-827-3631 , 913-864-4836 
SUBMODULE TITLE: Relative Gain and Antenna Pointing Angle Corrections 

SUBMODULE FUNCTIONS: This submodule Is designed to provide an output to be 

used to verify and locate possible errors in the relative gain table (Module 
1,6) and in the antenna pointing angle estimates. A maximum likelihood pro- 
cess will be applied using the means and standard deviations determined in 
Submodule 1.12.3 to estimate values for the antenna gain bias and pointing 
angle for each combination of beam, cell, and polarization. Average biases 
and pointing angle for each beam and polarization are also computed. 

INPUTS: 

(l) Mean a 0 , standard deviation of a 0 , and average Incidence angle for the 
Amazon rain forest passes (Submodule 1.12.3 outputs), 

OUTPUTS: 


Key Parameters No. of parameters 

Antenna biases (a) 600 

A 

Pointing angles (9 p ) 600 

Average biases (a) 12 

Average Pointing angles ( 0 p ) 12 

AUXILIARY DATA: Nominal Pointing Angle, (0^ ) Aa, A0 

EXCEPTIONS: TBD ’ P 

PROCESSING REQUIRED: 

For each combination of beam (i), cell (j), and polarization (k) perform the 
following algorithm (i x j x k = 600 times). 


(a) a Q = 1.0, 0 Q p = 0p (the nominal pointing angle) 

(b) Determine a 3 x 3 matrix, g as follows: 

m,n 

°m = °o + m Aa; m = , ’ 0 ’ 1 
9 n p = 0 o p + n A0p; n = 


-1 

9 m,n ~ 2 J, 


{a° _ a [G/Gg( ; e n p)] i2 

i,j,k,2 m [G/G0 (0 )]2 A 

P 


- a "^rr " 

m r ^ ~ i n U A J 
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M2.7(2) 


where 

G/G0 ( ) Is the output of Submodule 1 .12,5 at *0. , . 0 
o° ( ) Is the output of Submodule 1.12,6 at "©j j k ^ 
a 0 and O’ are the output of Submodule 1.12.3 

(c) Search the matrix g m>n for the maximum value 9 mmax ^ max 
If mmax * nmax “ 0 go to (d); 

otherwise % 

S o p ' 0p nmax 

repeat (b) 

(d) Compute the constants of a bivariate Interpolating polynomial for g. 


1 .0 4. 

a ■ ~z~ ■ s o,o + ■ 

b = iL0 + lU0 
2 2 

g 0 -i 

C = —I 9 0 f 0 + 

- ~ 9 °.-i + 9 °.' 


9 1,0 

2 


9 0,1 

2 


9 0,0 " 9 1,0 " 9 0,1 + 9 1,1 



where g are the points of the 3x3 matrix 
m , n 


(e) Solve for antenna bias* pointing angle 


50 a 

J»1 p JJ»k 


1.12.7(3) 
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XI. COMMENTS; 

(1) It Is assumed that this processing must take place only for the outer 
50 cells. 

(2) Initial values for Act ,A0p are 0.2 and 1.0 respectively, 

(3) A test should be made to tell If there are at least « 10 elements In 

each beam, cell, polarization bln of the input data before this submodule 

Is called. The bins are set to handle up to 20 elements to allow for 

the varying rates of filling. (On the long term average the bins will 

fill at the same rate). After the first 10 levels are processed by this 

submodule, further data from 1.12,3 may be stored over this data, This 

test can be performed by Interrogating P, , of Submodule 1.12,3. 
o o ,J 

(*♦) Mean a and standard deviation of a should be In ratio form. * 
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1 . 12 . 8 ( 1 ) 


I. SUBMODULE REFERENCE NO: 1,12.8 
II. SENSOR: SCATT 

III, ORIGINATOR/PHONE: E.M. Bracalente, I.J. BJ rrer/804-827-3631 , 913-864-4836 

IV. SUBMODULE TITLE: Long-Term Monitoring of Relative Biases 

V. SUBMODULE FUNCTION: Given the value of the pointing angle for each beam and 

polarization, compute the relative bias for each of the cells In that beam; 
also, provide long term checks for both biases and antenna patterns. 

VI. INPUTS: 

(1) Mean NRCS (o°) and average Incidence angles for the Amazon rain 
forest passes (SUBMODULE 1,12.3 outputs) 

(2) Previously determined average pointing angles (0 p ) for each beam (l) and 
polarization (k) (Submodule 1.12.7 outputs). 

(3) Initial guess for relative biases - outputs from either submodule 
1.12.7 («) or some first guess (e.g., 1,0). 

VII. OUTPUTS 

Key Parameter No. of Parameters 

Antenna bias (a) 600 

Average biases s(cx ) 12 


VIII. AUXILIARY DATA: Standard radar target, or output from Submodule 1.12.4, Aot 

IX. EXCEPTIONS: TBD 

X. PROCESSING REQUIRED: 


4 


* 


For each combination of beam (l), cell (j) and polarization (k) 
perform the following algorithm (I x j x k “ 600 times). 

(a) a = a, , If available; otherwise a =1.0- 

O I y J f K O 


(b) 


Determine a 3 element vector, g , as follows: 

’ m 

« m = a„ + mAa; m «* 1,0,1 
mo ’ 


’m 


1 

2 


10 

I 

£=1 




U,M 


" Va ]) ‘ 


where a°( ) Is the output of submodule 1.12.6 at ®j j ^ £ 
o° and 0 are the output of sumbodule 1.12.3 
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M2. 8(2) 


A , 9i ” g w i 

(c) a, , ,, - cu ** Aa(- - — j— r-) 

«»J» k 0 - 2g Q + gj 

1 50 * 

(d) «| f k * W jJj a l ,J ,k 
XI. COMMENTS s 

(1) It Is assumed that this processing will take place only for the 
outer 50 cells of the Scatt Beams. 

(2) Initial valM? for Aa Is 0,2 

(3) As In Submodule 1.12.7* a test must be made to be sure that there 
are 10 elements In each bln. (For further details see Submodule 
1.12.7 see XI 3). 

(h) Mean 0° and standard deviation of o° must be In ratio form. 


